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Abstract 



tic limit and derive the limit stochastic Schrodinger equations for the wave 
function conditioned on indirect observations using only the von Neumann 
projection postulate. We show that the diffusion (Gaussian) situation is 
.^j. , universal as a result of the central limit theorem with the quantum jump 

f^ ■ (Poissonian) situation being an exceptional case. It is shown that, start- 

ing from the correponding limiting open systems dynamics, the theory 
of quantum filtering leads to the same equations, therefore establishing 
consistency of the quantum stochastic approach for limiting Markovian 
models. 



qh' 1 Introduction 



The problem of describing the evolution of a quantum system undergoing con- 
^^ , tinual measurement has been examined from a variety of different physical and 

Vh ' mathematical viewpoints however a consensus is that the generic forms of the 

stochastic Schrodinger equation (SSE) governing the state tp^ {uj), conditioned 
on recorded output ui, takes on of the two forms below: 

MVt) = {L-\t)\ij,)dqt+(-iH-^{L^L-2XtL + X't)')\^,)dt, (1) 

\d^t) = (^^=^) 1^*) '^"^ + {-'^ - \^^^ - ^^' + ^/^^) 1^*) ^*(2) 

Here H is the system's Hamiltonian and L a fixed operator of the system which 
is somehow involved with the coupling of the system to the apparatus. In (^ we 

have At [uo) = - (i/^j (w) | (L^ + L) ip^. (w)) and {qt ■ t > 0} is a Gaussian martin- 
gale process (in fact a Wiener process). In (jSJ, i^t (uj) — (ipt (uj) \ L^Lip^. (cj)) and 



{fit : t > 0} is a Poissonian martingale process. The former describing quantum 
diffusions [1-7], the latter quantum jumps [8-10]. 

(By the term martingale, we mean a bounded stochastic process whose cur- 
rent value agrees with the conditional expectation of any future value based on 
observations up to the present time. They are used mathematically to model 
noise and, in both cases above, they are to come from continually de-trending 
the observed output process.) 

There is a general impression that the continuous time SSEs require ad- 
ditional postulates beyond the standard formalism of quantum mechanics and 
the von Neumann projection postulate. We shall argue that this is not so. Our 
aim is derive the SSEs above as continuous limits of a straightforward quan- 
tum dynamics with discrete time measurements. The model we look at is a 
generalization of one studied by Kist et al. ^21 where the environment consists 
of two-level atoms which sequentially interact with the system and are subse- 
quently monitored. The generalization occurs in considering the most general 
form of the coupling of the two level atoms to the system that will lead to 
a well defined Markovian limit dynamics. The procedure for conditioning the 
quantum state, based on discrete measurements is given by von Neumann's 
projection postulate. Recording a value of an observable with corresponding 
eigenspace-projection 11 will result in the change of vector state -0 i-^ p^^^^Ilip 
where p = (?/'|n-0) is assumed non-zero. Let us suppose that at discrete times 
t = T, 2t, St, . . . the system comes in contact with an apparatus and that an in- 
direct measurement is made. Based on the measurement output, which must be 
viewed as random, we get a time series [ipj) . of system vector states. The ques- 
tion is then whether such a time series might converge in the continuous time 
limit (r — *■ 0) and whether it will lead to the standard stochastic Schrodinger 
equations. We apply a procedure pioneered by Smolianov and Truman jJSI to 
obtain the limit SSEs for the various choices of monitored variable: a key feature 
of this procedure approach is that only standard quantum mechanics and the 
projection postulate are needed! 

The second point of the analysis is to square our results up with the the- 
ory of continuous-time quantum filtering 01,21], m],2S|. This applies to uni- 
tary, Markovian evolutions of quantum open systems (that is, joint system and 
Markovian environment) described by quantum stochastic methods J17 | .jl8 | .jl9 | . 
Our model was specifically chosen to lead to a Markovian limit. Here we show 
that filtering theory applied to the limit dynamics leads to exactly the same 
limit SSEs we derive earlier. Needless to say, the standard forms (QJ and ((SJ 
occur as generic forms. 

We show that if the two level atoms are prepared in their ground states then 
we obtain jump equations (0) whenever we try to monitor if the post-interaction 
atom is still in its ground state. In all other cases we are lead to a diffusion 
equation which we show to universally have the form ^ . 



2 Limit of Continuous Measurements 

Models of the type we consider here have been treated in the continuous time 
hmit by (lOlimi ^^'^ E2I- ^^ ^^^^ section we recaU the notations and results 
of |22| detailing how a discrete-time repeated interaction-measurement can, in 
the continuous time limit, be described as an open quantum dynamics driven 
by quantum Wiener and Poisson Processes. 

Let Hs be has state space of our system and at times t — t, 2t, 3t, ... it 
interacts with an apparatus. We denote by Ti.E.k the state space describing the 
apparatus used at time t = kr - this will be a copy of a fixed Hilbert space Ti.E- 
We are interested in the Hilbert spaces 

[t/r] oo 

^*i^^nEM, ^f^'= (g) nE,k, ^(^)=.Ff^®4^^ (3) 

k=l k=lt/T\+r 

where \_x\ means the integer part of x. (We fix a vector bq G He and use this 
to stabilize the infinite direct product.) We shall refer to J^^P and JF.j as the 
past and future environment spaces respectively. 

We are interested only in the evolution between the discrete times t = 
T, 2r, 3t, . . . and to this end we require, for each fc > 0, a unitary (Floquet) 
operator, Vfc, to be applied at time t = kr: its action will be on the joint space 
Hs (X" ^'-'^^ but it will have non-trivial action only on the factors Hs and HE,k- 
The 14 's will be copies of a fixed unitary V acting on the representative space 

(t) 

Hs ® He ■ The unitary operator U^ describing the evolution from initial time 
zero to time t is then 

ui^^^Vit/rr--V2Vi (4) 

It acts on Hs'^J^^] but, of course, has trivial action on the future environment 
space. The same is true of the discrete time dynamical evolution of observables 
X £ ^ {Hs), the space of bounded operators on Hs, given by 

4^Hx) = U^^^Hx^lr)ut\ (5) 

where Ir is the identity on J^'"^^ . The discrete time evolution satisfies the dif- 
ference equation 



[/l^/Vj +, - C/l^/Vj j = (nt/rj +. - 1) ^Lt/Vj • (6) 

The state for the environment is chosen to be the vector ^'^'^^ on J^^'^'> given by 

**^' = eo (g) eo (g) eo (g) eo • • • 

and, since cq will typically be identified as the ground state on He, we shall call 
^^'^> the vacuum vector for the environment. 



2.1 Spin- 1 Apparatus 



2 

For simplicity, we take He = C^. We may think of the apparatus as comprising 
of a two-level atom (qubit) with ground state eo and excited state ei. We 
introduce the transition operators 

a+ = |ei) (eol cr^ = |eo) (ei| 

The copies of these operators for the fc— th atom will be denoted a^ and a'^ . The 
operators aj^ are Fermionic variables and satisfy the anti-commutation relation 

{a±,a±} = 0, KT,C7+} = 1 (7) 

while commuting for different atoms. 

2.2 Collective Operators 

We define the collective operators A^ {t; r) , A {t; r) to be 

Lt/rJ [t/rl 

A±(i;r):=V7 5^a±; A (t; r) := J] a+a". (8) 

fc=i fc=i 

For times t, s > 0, we have the commutation relations 

'tAs 



[A-{t;T),A+{s;T)] 



T 



2rA(t As,t), 



[A (i; t) , y4± (s; t)] = ±A± (i A s; r) 



where s At denotes the minimum of s and t. In the limit where r goes to zero 
while s and t are held fixed, we have the approximation 

[A- {t;T),A+{s;T)] w i A s. (9) 

The collective fields A^ (i; t) converge to a Bosonic quantum Wiener processes 
A^ as r ^- 0, while A (i; r) converges to the Bosonic conservation process At 
|17|. This is an example of a general class of well-known quantum stochastic 
limits 122], mi- Intuitively, we may use the following rule of thumb for t = jr : 

T~dt, ^aJ~dAt, 

V7a+ ~ dA+, a+aj ~ dA*. ^"""^ 

These replacements are usually correct when we go from a finite difference equa- 
tion involving the discrete spins to a quantum stochastic differential equation 
involving differential processes. 

The limit processes are denoted as Aj^ — A^ , A'}^ = A^ and Al^ = At 
respectively and we emphasize that they are not considered to live on the same 
Hilbert space as the discrete system but on a Bose Fock space r+ (i^ (R^, dt)) . 
(See the appendix for details.) We also set ylj° = t. 



2.3 The Interaction 

The fc— th Floquct operator is then taken to be 

where 



H, 



(r) 



1 , _ 1 ,1 

-Hii (g) (tIcj^. + —^HiQ ^(71 + —i=Hqi (g) (jj^ 



Hn 



(11) 



(12) 



Here we must take Hn and Hqo to be self-adjoint and require that (ii^oi) — Hiq. 
We may identify Hqq with the free system Hamihonian Hs- We shaU assume 
that the operators Haf3 are bounded on Tis with Hn also bounded away from 
zero. 

The scaling of the spins a^ by t^^/^ is necessary if we want to obtain a 
quantum diffusion associated with the Hiq and Hqi terms and a zero-intensity 
Poisson process associated with Hn in the t —^ limit. 

We shall also employ the following summation convention: whenever a re- 
peated raised and lowered Greek index appears we sum the index over the values 
zero and one. With this convention, 



Hi^^ ^ Hep 





u. 





(13) 



were we interpret the raised index as a power: that is, [x] =1, [x] = x. 



2.4 Continuous Time Limit Dynamics 

We consider the Floquet unitary on Jis ® "He given by 



-1/3' 



V 



exp < —irHaf} ® 




ex 


a 


1 + tL^p (g) 




a 


O 


7?_ 
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where ~ means that we drop terms that do not contribute in our r ^ limit. 
Here the so-called Ito coefficients L^p are given by 



L. 



a/3 



-iH 



that is, 



ill 



-iHi 



aP 



L 



7 , —Hal (Hn)" Hip, 



(14) 



n>2 



-iHi 



10 



-iHi 



Lqi — Hoi- 



Hii 



Lo 



Hu 
-iHoo + Hoi 



-Hio] 



-iHi 



1 - 1 + JHu 
{Hiif 



H 



10- 



The relationship between the Hamiltonian coefficients Hajs and the Ito coeffi- 
cients was first given in j53] . We remark that they satisfy the relations 

L<.p + 4. + 4.^1/3 = 0. (15) 

guaranteeing unitarity '17' . Consequently we have 

Ln^W-l; Lio = L; Loi = -L^W; Loo = -iH - ^L^ L (16) 



with W — exp {—iHii} unitary, H — i^oo^-ffoi^^Tfp^T^^io self-adjoint and 



-H"ii-sin(_ffii) 

L is bounded but otherwise arbitrary. (Note that ^ ^^^ > for a; > 0.) 

2.4.1 Convergence of Unitary Evolution 

In the above notations, the discrete time family < U^ '■ t> 0> converges to 

quantum stochastic process {Ut : i > 0} on Hs <Xi r+ (i^ (M+, dt)) in the sense 
that, for all u,v G Tis, integers n,m and for all (pjjipj G L^ (Es^,dt) Riemann 
integrable, we have the uniform convergence (r -^ 0) 

(a+ (</>„„ r) •• •A+ (01, r) w® vI;M| f/M A+ (^„,t) •• •A+ (^i,r) z;0 vI^M^ 

-. {A+ {c^J ■■■A+ (0i) u<»^Ut A+ (V'„) ■■■A+ (V-i) V (g> *) (17) 

The process Ut is moreover unitary, adapted and satisfies the quantum stochastic 
differential equation (QSDE, see appendix) 

dUt = L^p ® dAf Ut, Uo = 1. (18) 

2.4.2 Convergence of Heisenberg Dynamics 

Likewise, for X a bounded observable on Tis the discrete time family I Jj {X) > 

converges to the quantum stochastic process Jt {X) = UJ {X (g) 1) Ut on Hs 'Si 
r+ (L^ (K^, dt)) in the same weak sense as in ()17|l . The limit Heisenberg equa- 
tions of motion are 

dJtiX) = JtiC^pX)sdAf, JoiX)=Xsl (19) 

where 

Cc^p (X) := lI^X + XL^p + L\^XLip. (20) 

We remark that Cap (1) = 0, by the unitarity conditions (I15II . A completely pos- 
itive semigroup {St : i > 0} is then defined by {u\ Et {X) v) := (u (g) ^| Jt {X) w (g) *) 
and we have St — exp{i£oo} where the Lindblad generator is Loo {X) = 
\ {L\X\ L + iit [x,L] - i [X,H] with L = "'^"^[[^^ Hw and H = Hoo - 

rj -ffii— sin(ffii) TT 

^01 — (7?r7? — ^10- 

We remark that such QSDEs occur naturally in Markovian limits for quan- 
tum field reservoirs [21], [2H1j see also |I^ . 



3 Conditioning on Measurements 



We now consider how the measurement of an apparatus indirectly leads to a 
conditioning of the system state. For clarity we investigate the situation of a 
single apparatus to begin with. Initially the apparatus is prepared in state eo 
and after a time r we have the evolution 



(f) ig) eo ^ V {4> ig) Co) c^ (l + tLqo) (g) eo + ^/rLiocf) (g) ei. 
We shall measure the o-^-observablc. This can be written as 



(21) 



a 



a 



a 



|e+; ;e+| 



e- 



where |e+) = ^ l'^i) + 7f l^o) and |e_) = ^ l'^i)~7f l^o)- (Actually, the main 
result of this section will remain true provided we measure an observable with 
eigenstates |e±) different to |eo) and |ei). We will show this universality later.) 
Taking the initial joint state to be (g Cq, we find that the probabilities to 
get the eigenvalues ±1 are 



p± = (0 (8) eol F^ (Is «> n±) F </) (g) eo) 
where we introduce the expectation 



V2 



[l ± 2AV7] +o[t 



Jil2 



A = i (</)| (L + Lt) 0) . 

A pair of linear maps, V± on ()s arc defined by 

(V±(/)) (g) e± = (Is ® n±) F (0 (g eo) 
and to leading order we have 



(22) 



^-^T. 



i±LV7- 



1 



-iH- -UL 



The wave function W±0, conditioned on a ±-measurement, is therefore 



(23) 



>V± will be non-linear as the probabilities p± clearly depend on the state (/). We 
then have the development 



W±</)~0±\/7(i-A)0 + r 



1 l'\ 
-iH- -L^L + \ -A-L 

2 V2 



(24) 



We now introduce a random variable rj which takes the two possible values 
±1 with probabilities p±. We call rj the discrete output variable. Then 



EM = P++P- =2AV7+0(t) 
E[?7^] = p++P- = 1 + 0{t). 



(25) 
(26) 



Now let us deal with repeated measurements. We shall record an output 
sequence of ±1 and we write ry. as the random variable describing the jth output. 
Statistically, the 77 ■ are not independent: each ry will depend on r]^,- ■ ■ , rjj-i- 

We set, for j = [t/T\ and fixed (f) G ^s, 



/.("'= v^„---v^^0, vr^=w^^...w^^0 



1 



We shaU have Pr {rj^^^ = ±1 } ~ i 1 ± V^A 
and 



(r) 
J 



, where A:- = 



/,(^) 



(r) _ 



bV- 



(27) 



i(^(^)|(L + Lt)V.M), 



Vlo+i\ 



2A(")^/7- 



0(r) 



EI 



{v,+if = l + 0(r) 



(28) 
(29) 

rjjY The state 

V'(I+i)i- after the (j + l)-st measurement depends on the state Vjr ^^^ '^j+i 
and we have, to relevant order, the finite difference equation 



where E^ is conditional expectation wrt. the variables (?7]^, • 






We wish to consider the process 



-iH 



> 



,(-) 



-A 



(r) 



L 



(r) 



^r 



g(^) it) = 



Lt/rJ 



however, it has a non-zero expectation and so is not suitable as a noise term. 
Instead, we introduce new random variables C- '■= Vj ^ V^'^^Z-i and con- 
sider q'-^^ (t) — V^^j=i Cj- We now use a simple trick to show that it is 
mean-zero to required order. First of all, observe that Wj_i [CA = and so 
E [CJ] - E [EJ_i [CJ]] - O (t). Similarly E [C|] =1 + {^). It follows that 
{q {t) : t > 0} converges in distribution to a mean-zero martingale process, 
which we denote as {qt '■ t> 0}, with correlation E [qtqs] = tAs. We may there- 
fore take qt to be a Wiener process. Likewise, {g^"^-* {t) '■ t > O} should converge 
to a stochastic process {qt '■ t > 0} which is adapted wrt. q: in other words, qt 
should be determined as a function of the q-process for times s < t for each time 
i > 0. In particular, 



qt = qt + 2 / Xsds 



where At 



{ijt\{L + L^)^Pt)- 



In the limit t — > we obtain the limit stochastic differential 



\di;t)^iL-Xi)\^Pt)dqt + 



-iH 



-L^L + A, I -At - i 



l^t) dt 



with initial condition |V-'o) = \4')- Ii^ terms of the Wiener process q we have 



\di^^) = {L-Xt)\i^t)dqt + 



iH ~-{L^L-2\tL + X 



l^t) dt. (30) 



This is, of course, the standard form of the diffusive Stochastic Schrodinger 
equation (Q. 



3.1 Universality of Gaussian SSE 

We consider measurements on an observable of the form 

X ^x+ \e+) (e 



- x^ e_ e_ 



(31) 



where x^ and x_ are real eigenvalues, while e^ and e_ are any orthonormal 
eigenvectors in He not the same as eo and ei. Generally speaking we will have 



e+ = y/qeo + e^^y^l - qei, e_ = v^l - qeo - V^^'^ei 



(32) 



with < q < 1. For convenience we set q+ = q and g_ = 1 — q. The phase 
9 G [0, 27r) will actually play no role in what follows and can always be removed 
by the "gauge" transformations eo ^^ eo, e ei ^^ ei which is trivial from our 
point of view since it leaves the ground state invariant. We therefore set 9 = Q. 
The probability that we measure X to be x± after the interaction will be 



^ '?±[l + ?7±2AV7-//(l-ry^)r] 
where we introduce the expectations 

v^{<ly\L^L4>) = \\L4' 
and the weighting 

V- = 



(33) 
(34) 



"-^^T,^ 



(35) 



We may now introduce a random variable 77 taking the values rj^. with proba- 
bilities p± . We remark that 



EW = p+v++P-V-=^^V^+0{t), 



[V^] = P+vl + P-V- = '^ - 2X- 



ql 



Vq+T- 



't + 0{t). 



(36) 
(37) 



We then have the finite difference equation 

r 1 X 

+T 



Jr)2 



-^-2^'^+2(i-^a^)J'^(;)+^a^)i^^. 



,M 



J^)2 



^Af^-A^i^ 



4'^ (38) 



which is the same as before except for the new term involving i/Z! = ( ip^^ \ L^ L t/j 



',Mi 



/,(^) 



We may replace the ry^-terms by their averaged value of unity when transferring 
to the limit r ^ 0: in particular the term with v^ is negligible. Defining the 



,(t) 



processes q^' and g^ ' as before, we are lead to the same SSE as (|5n|) . 



3.2 Poissonian Noise 

Now let us consider what happens if we take the input observable to be cr+cr^. 
(This corresponds to the choice q+ = 1, q^ — 0.) We now record a result of 
either zero or one with probabilities p^ — {(f)® e[j\V^ (1 (8) Ilg) V<t) ® cq) where 
e = 0, 1 and li^ — \e^) {e^\. Explicitly we have 



Pq 



1 - VT, pi 



As before, we define a conditional operator Ve on Hs by {Ve(f>)<E)ee = {I ® Us) V {(fx® eo) 

I/O 

and a conditional map We = {Pe) ^e- Here we will have 



Wo./): 



1 + T ( -iH- -L^L+-v 

2 2 



Wi(t>~-^L(t>. 



Iterating this in a repeated measurement strategy, we record an output series 
■' ^^'^' ■ ■ ■ ' of zeroes and ones with difference equation for the conditioned 



,£ 



2 ' 



states given by 






tP 



(r) 



Jr) 
'J + 1 



L 



^ 



ir) 



t1 -e 



J + 1 



-iH--L^L+-J''^ 
2 2 ^ 



V' 



{r) 



(39) 



where vp := ( ■4)p\L^ Lijjp \ . The ep 's may be viewed as dependent Bernoulli 
variables. In particular let Ej [•] denote conditional expectation with respect to 
the first j of these variables, then E^ exp 
now define a stochastic process ni'' by 



(iue^P^ ~ 1 - Tv'^p (ie" - 1). We 



Lt/rJ 

h — 2^ tj 



and in the limit r ^ this converges to a non-homogeneous compound Poisson 
process {rit : i > 0}. Specifically, if / is a smooth test function, then 



limE 



cxp ■ 



Vt/T\ 

E 



njr)e 



(r) 



= E 



exp 



ds 



Vs f e 



i/(s) 



10 



with mean square limit i^t '■= liiHr^o '^i f/ri • ^^^ Ito table for rit is (drit) — dnt, 
and we have E^j [dnt] = vtdt where E(] [•] is conditional expectation with respect 
to nj. 

The limit stochastic Schrodinger equation is then 



l#* 



L 



J^'^ \^t) dnt + (-iH - i {L^L - ut) ] |^,) dt (40) 



and one readily shows that the normalization ||'0t|| — 1 is preserved. Replacing 
Tit by the compensated Poisson process ht = rit ~ Jq Vsds, we find 



\d^t 



IV't) diit 



-iH-]^{L^L + vt)+vtL\\^t) dt (41) 



This is the standard jump- type SSE |(2J). 



3.3 Discrete Input / Output Processes 

In quantum theory the probabilistic notion of events is replaced by orthogonal 
projections. For the measurements of ax: the relevant projectors are 11^? = 



l±cr 



U) 



and so far we have worked in the Schrodinger representation. We 



note the property that, for j ^ k, 



n2\i4 



= 0. 



In the Heisenberg picture we are interested alternatively in the projectors 



n2^=t/|;)^n2^c/r) 



JT 



(42) 



Note that 



n2"\ni^) 



for j < k. 

The family <^ 11^ : j — 1,2,- ■ ■ > is auto-commuting: a property that is some- 
times referred to as leading to a consistent history of measurement outputs. 
The family < W^' : j — 1,2,- ■ ■ > is likewise also auto-commuting. To see this, 
we note that for n > j, 

^ nr ± ^ nr 



v^---v,lii'-£v„---Vi 
v^ ■■■ vj u^£ Vj---Vi--. 



n 



± 



nl\nl^ 



^ nr 



and so, for any j and fc, 

need only take n greater than both j and k. 

For a given random output sequence tj = (r^j^, 772, 



n2\n(^^' 



Unr — where we 
we have for n — [i/rj 



(^i^^ {ri)) 



.$ 



;) = (n('^)K)---(n(iVi)0®$(- 
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and the right hand side can be written ahcrnatively as 

n(") • • • n(i) c//^^ (j, ® $(^) or ul^^ M") • • • M^^ </) (g) $(^^ 

The probability of a given output sequence (r/^^, • • • , r/^) is then 



4^'> (r,) 



= {Xl 



J^) 






"*"^Ti "*" '?! -^^ 



n^^ • • • ul^^^ (j) ® $(") 



where xi^' := t/^^^V <» *^^^- 



Likewise, we have V't (j?) = 0t (''?) 
of the system we have the random expectation 



a(^) 



(?7) and for any observable G 



^ 



(r) 



G^r 



(r)\ _ 



dr) 



i>V 



Xi 



(r) 



G®i^^Aui-^---ni'} x^^ 



4) ® $(") jf ^ (G) ii(,'^) • • • ii(,i) (t> ® $(") 



Let us introduce new spin variables ct defined by 



~^t 



U^^^^ (a±) U^^^ 



In particular, let '^^ ~ '^t + S-- then tl^^' = ^ [l ± a^] . We may the construct 
the following adapted processes 



[t/r] 



A(t;r):= 



[t/r] 



The current situation can be described as follows. Either we work in the 
Schrodinger picture, in which case we are dealing in the quantum stochastic 
process Q {t; r) = A~^ {t; r) + A^ {t; t), or in the Heisenberg picture, in which 
case we are dealing with Q {t; r) = U^^^'' X (i; r) U^''^ = A+ {t; t) + i" (i; r). 
Adopting the terminology due to Gardiner, the Q-process is called the input 
process while the Q is called the output process. We may likewise refer to the 
n^ s as input events and the n^i s as output events. 

It is useful to know that, to relevant order, the output spin variables are 

= 4-i)r(^)^''J + 4-i)r(^)^ (43) 

We shall study the continuous-time versions of these processes next. 
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4 The Stochastic Schrodinger Equation 

We now review the simple theory of quantum fihcring. Let Y — {Yt : t > 0} be 
an adapted, self-adjoint quantum stochastic process having trivial action on the 
system. In particular we suppose that it arises as a quantum stochastic integral 



Yt=yo+y2 f Ya,f3 (t) dA 



Q/3 





where the Yap (t) = {Yf^a (t)) are again adapted processes and yo G C. We 
shall assume that the process is self-commuting: 

[Yt,Ys] =0, for alH,s > 0. 

This requires the consistency conditions [Yajs (t) ,Ys] — whenever t > s. 

The process Y can then be represented as a classical stochastic process 
{yt ■ t > 0} with canonical (that is to say, minimal) probability space (ri,S,Q) 
and associated filtration {Sj] : t > O} of sigma-algebras. For each i > 0, we 
define the Dyson-ordered exponential 

f cxp \ f f{u)dYA:=y2- f dYt^ ■ ■ ■ dYt, f (i„) • • • / (t,) 



n>0 



where A„ (i) is the ordered simplex {(i„, • • • , ti) : i > t„ > • • • > ii > 0}. The 
algebra generated by such Dyson-ordered exponentials will be denotes as £^. 
Essentially, this is the spectral algebra of process up to time t and it can be 
understood (at least when the y's are bounded) as the von Neumann algebra 
C^ = {Ys : < s < t} where we take the commutants in *B (Hs <Xi ^t])- In 
the following we shall assume that €q^ is a maximal commuting von Neumann 

sub-algebra of *B {J^t])- In other words, there arc no effects generated by the 
environmental noise other than those that can be accounted for by the observed 
process. The commutant of (£X will be denoted as 

2[^ = (^t^)' = {^e*B(Hs®J^) : [z,A] =o,vze e:J[} 



and this is often referred to as the algebra of observables that are not demolished 

-t] 



by the observed process up to time t. We note the isotonic property €Yi C 2^!^ 



whenever t < s and it is natural to introduce the inductive limit algebra £^ 

limt_^o£^- 

From our assumption of maximality, we have that 

2t^ = » (Hs) ® C^ ® *B (.F(t) . 



A less simple theory would allow for effects of unobserved noises and one would 
"X^ as the c( 
l^^intoC^, 



have £^ as the centre of 'QiK . One is then interested in conditional expectations 
from 2l?[ into €X . Here we are interested in the Hilbert space aspects and so we 
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take advantage of the simple setup that arises when C^ is assumed maximal. 
(For the more general case where €j, is not maximal, see IIS).) 

It is convenient to introduce the Hilbert spaces TiX and G^ defined though 



21^ {Hs ® ^t]) ^ n\ ® {C*(,} , €l (vft]) ^ Gl 

where we understand TiX as a subspace of Tis ® Tt\ and G\ as a subspace of 
T^, . From our maximality condition we shall have 

n\ = Us®G\. 

(Otherwise Jis ® G\ would be only a subset of H|^.) A Hilbert space isomor- 
phism 3( from l-is ® G\ to Hs ® L^ (f^, S, Q) is then defined by linear extension 
of the map 

(/. ® f exp j /" ](u) dvA ^ 0f exp | /" f (u) dy, 

where the Dyson-ordered exponential on the right hand side has the same mean- 
ing as for its operator- valued counterpart. 

We remark that, in particular, we have the following isomorphism between 
commutative von Neumann algebras: 

By extension, €.^ can be understood as being isomorphic to L°° (fi, E, Q). Now 
fix a unit vector (j) in Tis and consider the evolved vector 

which will lie in Tis (^GJt,- In particular, we have a E-measurable function 0j (•) 

corresponding to "^tXt- Here (|)^ is a Hs-valued random variable on (fi, I],Q) 

which is adapted to the filtration {Ej] : t > O}. We shall have the normalization 

condition 

" \cl>^[u)\\l'[l[du;] = l. 
Jn 

In general, ||(/)j (ti-')!!^ is not unity, however, as it is positive and normalized, we 
may introduce a second measure P on (17, E) defined by 

V[A]:^ f U,iu;)\\lQ[dco] 

J A 

whenever A dTjt. We remark that, for B E €X, we have 



iXtlBXt)^ htB3i^V[duo] 
Jn 
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It is convenient to introduce a normalized Ti^-valued variable V't defined 
almost everywhere by 



^ti^) := Uti^)\\s^ <Pti^) 



We now define a conditional expectation £X from 2t?[ to the von Neumann 



-t] ii^^i^t] 

l-n r -i- In /-I 4-<-~i 1 I <-~m t i n r-c i /-I <-~i n 4- i-fn .^ i-i + 1 <-\ vi 



sub-algebra €X by the following identification almost everywhere 



This expectation leaves the state determined by Xt invariant: 

Xt\£l][A]xt) = Jji^A^tA3^'^,)F[cLo] 

= iXtlAxt)- 

This property then uniquely fixes the conditional expectation. If we consider 
the action of £^ from £^, only, to £^ then this must play the role of a classi- 
cal conditional expectation E^, from S-measurable to Sj-measurable functions, 
again uniquely determined by the fact that it leaves a probability measure, in 
this case P, invariant. We have the usual property that E^, o E^, = E^^^, . We 
shall denote by E^ = E^, the expectation wrt. P. 

Let us first remark that the classical process {yt : t > 0} introduced above 
is not necessarily a martingale on (f7,S,P) wrt. the filtration |Sf] : t > O}. 
Indeed, we have 

Ey[dyt] = (XtldYtXt) 



where we define the output process Y by 



Yt := Ul Yt Ut 



From the quantum Ito calculus, we obtain 
rt jv^. TT I ^rrt 



dYt = WtdYtUt+dUjYtUt + UjYtdUt+dUjYtdUt 
+dUl dYt Ut + UJ dYt dUt + du} dYt dUt 
= U} Y^ti (t) Ut dAf + Ul U^ (1) Yt Ut dAf 

+U} (l\^Yip + Ya^iLip + iLYiiii^) Ut dAf. 



Noting that Cap (1) = 0, we see that 
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dYt = U} (y^p (i) + L\^Yip (y) + r„i (t) Lip + lLFh (t) ii^) C/^ dA,"^. 

In particular, we define A^ := UjA"Ut and they are explicitly 

dAt = di," = dAt + Jt {W^L) dA+ + Jt {L^W) dA^ + Jt [L^ L) dt 
dAt = dA\^ = Jt {W^) dA+ + Jt (L^) dt 
dA^ = dAf = Jt [W) dA- + Jt {L) dt 

with i™ = t. 

We remark that & [dyt] = ytdt where 

yt = {xt\ (>^oo (i) + ilo^io (t) + >oi (t) Lio + l1o>1i (^) ^lo) Xt) 

and we use the notations yap (t; •) = 3ti^a/3 (t) ^3^ ■ Therefore, a martingale on 
(17, 1],P) wrt. the filtration {S^] : t > O} is given by the process {ijt ■ t > 0} 
defined as 

dyt i^) = dyt (lu) ~ {^t (^) I [-^^]" [-^l'^ V"* (^)) Vafi {t;uj) dt. 

4.1 Filtering based on observations of Qt = Af + A^ 

Let us choose, for our monitored observables, the process Qt = Af + A'[ . Here 
the output process will be Qt with differentials 

dQt = Jt {W^) dA+ + Jt (W) dA^ + Jt {L^ + L) dt. 

By the previous arguments, we construct a classical process y = q giving the 
distribution of Q in the vacuum state: as is well-known, this is a Wiener process 
and (SI, S, Q) will be the canonical Wiener space. (In fact, Jt is then the Wiener- 
Ito-Segal isomorphism 18 .) The corresponding martingale process will then be 
q defined through 

dqt = dqt - 2Xtdt, qo = 

where 

A, (cj) := i (^, {iv) I {L + Lt) v., (c.)) . 

A differential equation for ipt can be obtained as follows. The state Xt—^t <t>® 
^ will satisfy the vector-process QSDE 

dxt = LapdAf Xt = LaodAf Xt (44) 

since we have dAf^Xt — UtdAf^(j)iS) ^ = - that is the Ito differentials dAt and 
dA^ commute with Ut annihilate the Fock vacuum. It is convenient to restore 
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the annihilation differential, this time as LiqcIA^ (f>^ = 0, in which case we obtain 
the equivalent QSDE 

dxt = -Uh+ -L^Lj Xtdt + L dQt Xf 

It should be immediately obvious that the process </)((•) will satisfy the sde 
\d(j)f.) — L \(f>f) dqt — (iH + \L^ L) \(f>^) dt. Here which we shall write (j)^ (•) as 
\<|)^ (•)) to emphasize the fact that it in Ti^-valued process. From the Ito rule 
{dqt) ~ dt, we find that 

The derivative rule is 

dUtV = iUtf+dUtfy'^-iUtf^'"^ 



^Erf)ii^*r"('^ii^*iif (45) 



"*" ^\ k 

fe>l 



and here we must use the Ito rule d \\4't\\ ~ 2Xtdt where here At := ^ (V'tl {L^ + L) tpt)- 
This leads to 

dUtir' = -lutir' >^tdqt + lut\\-' x'tdt- 

This leads to the SDE for |V't): MV't) = ll^tlP^ \d(t)t) + d (\\4>t\\^^ 
d [\\4>t\\ ) \d(j>t) and this is explicitly 



MVt) = (i - Xt) Ht) dqt +{-iH- h^L - XtL + ^A? ) I^J dt. (46) 



Finally, substituting in for the martingale process q we obtain 

\d^t) = {L- At) l^t) dqt + (-iH - i {L^L - 2\tL + A?)") I^J dt. 

4.2 Filtering based on observations of Aj 



(47) 



Let us now choose, for our monitored observables, the gauge process Af. Unfor- 
tunately, we hit on a snag: the gauge is trivially zero in the vacuum state, that 
is, it is a Poisson process of zero intensity. A trick to deal with this is to replace 
the gauge process with a unitarily equivalent process Af given by 

A{:=e^"WV^^(/)Ate^+(/)-^"(/) 

for / e L^ (R+, rft) a real- valued function with / (t) > for all times i > 0. The 
process is defined alternatively by 

dA{ = dAt + f (t) dA+ + f (t) dA-; + / {tf dt, aI = 0. 
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It satisfies the Ito rule dA{ dA{ = dA{ and we have that /^| dA{ ^\ ^ f {t)'^ dt. 

We see that A( corresponds to a classical process y — n^ which is a non- 
homogeneous Poisson process with intensity density /^ and we shall denote by 
(17, E, Q) the canonical probability space. 

Now, from (gH), we find dxt = - {iH + jL^ L + fL) Xt dt + f^^L dn{ Xt and 
the corresponding Tig-valued process satisfies 



1 



iH+-L^L + JL 



dt + f-^L 



K) dut 



f 



from which we find that 

dUtf = -^{^t\{L^ + f){L + f)^,) (dn{-f{tfdt) 

Substituting into H45(l we find after some re-summing 

f{t) 



d\\4>t 



^^yt + 2f{t)\t + f{t)' 



-l\ dn{ 



where vt (f^) '■— {i't {^) I ^^Ltp^ (w)) and At (uj) is as defined above. Note that 
ut + 2/ (i) Xt + f {tf = (^, [u] I (Lt + / (<)) {L + f (t)) Vt [to)). The resulting 
sde for the normalized state ip^ is then 



M,.) ^ l'-^I^^^±llMt±l^]\,,)dnf 



^ut + 2/ [t) Xt + f (tf 

' - ^L^L -f{t)L+^ [vt + 2/ (i) At]) IVt) dt. 



Now n^ is decomposed into martingale and deterministic part according to 

dn^ = dhf -f (ut + 2/ (i) At + / {tf\ dt 
and so we have 

L + !{t)- ^yt + 2f{t)Xt + !{tY 



\d^t 



iV-t) dn{ 



^Ut + 2f{t)Xt + J{tf 
iH - h^L -f{t)L+^[iyt + 2/ it) Xt] 

+ (L + fit)- ^^t + 2fit)Xt + fitf] ^{ut + 2f{t)Xt + f{tY 



l^t) dt. 
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We now take the limit / ^ to obtain the resuh we want and this leaves 
us with the sde 



^ J[' ) iV't) dnt + (-iH - h^L + ^i^t) \^,) dt. 



Here nt will be a non- homogeneous Poisson process with intensity i^t. 

5 Appendix 

5.1 Bosonic Noise 

Let 7i be a fixed Hilbert space. The n-particle Bose states take the basic form 
4>i'^ ■ ■ ■ ®4>n — X]o-e6„ '^cr(i) ® ■ ■ ■ ® 't'a-(n) whcrc wc sum over the permutation 
group 6„. The n-particle state space is denoted Ji^^ and the Bose Fock space, 
with one particle space H, is then r+ (H) := ©i^o^*^" '^^^^ vacuum space 
TC^'^ spanned by a single vector ^P. 

The Bosonic creator, annihilator and differential second quantization fields 
are, respectively, the following operators on Fock space 

^+ (V;) 0i(8) ••• «)(/)„ = Vn+ I tp(g)(l)i(g) ■ ■ ■ (g)<j>„ 

A~ (i/i) 0;^(g) ••• «)(/)„ = ^ ^(V'|(/')«)(/)i«)---(8)i^«) •••«)(/)„ 

dr (T) <^i0 • • • (^(/)„ - ^ 01® • • • ® (T0,) ® • • • ®0„ 

j 

where i^ eH and T e ^ (H). 

Now choose H ^ L'^ (R+, rfi) and on the Fock space T = r+ (L^ (R+, rft)) 
set 

A±:-A±(l[o,t]); A* := rfF (ip,*]) (48) 

where l[o,t] is the characteristic function for the interval [0,t] and l[o.t] is the 
operator on L^ (K^, dt) corresponding to multiplication by Iro.ti- 

An integral calculus can be built up around the processes A^ , At and t and 
is known as (Bosonic) quantum stochastic calculus. This allows us to consider 
quantum stochastic integrals of the type /„ {-Fio (t) <S) dAf + Foi (t) (g) dA^ + 
Fii (i) (g) dAt + Foo {t) (g) dt} on Ho «> r+ (L^ (M+, dt)) where Ho is some fixed 
Hilbert space (termed the initial space). 

We note the natural isomorphism r+ (L^ (K^, dt)) = JF^j (g JF(j where J^t] = 
r+ (i^ ([0,t] ,dt)) and T^t = T+ (i^ ((t, oo), di)). A family {Ft)^ of operators 
on Ha g) r+ (i^ (R^, dt)) is said to be adapted if Ft acts trivially on the future 
space H(t for each t. 
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The Leibniz rule however breaks down for this theory since products of 
stochastic integrals must be put to Wick order before they can be re-expressed 
again as stochastic integrals. The new situation is summarized by the quantum 
Ito rule d {FG) = (dF) G + F [dG) + {dF) (dG) and the quantum Ito table 



X 


dA+ 


dA 


dA' 


dt 


dA+ 














dK 


dA+ 


dA 








dA- 


dt 


dA- 








dt 















It is convenient to denote the four basic processes as follows: 



A't 



At, (1,1) 

^(3_} A+, (1,0) 

AT, (0,1) 

i, (0,0) 



The Ito table then simplifies to dA" dA^'^ = except for the cases 

dAfdAf = dAf. (49) 

The fundamental result '17' is that there exists an unique solution Ut to the 
quantum stochastic differential equation (QSDE) 

dUt = L^fi (g) dA^^, Uo - 1 

whenever the coefficients Lap are in *B(7io)- The solution is automatically 
adapted and, moreover, will be unitary provided that the coefficients take the 
form l(TH|l . 
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